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Demande internationale n° 

PCT/FR 98/00735 


Date du depot international (jour/mois/annee) 

10/04/1998 


(Date de priori te (la plus ancienne) 
(jour/mois/annee) 

14/04/1997 


Deposant 
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Le present rapport de recherche internationale, etabli par ('administration charges de la recherche internationale, est transmis au 
deposant conformement a I'article 18. Une copieen est transmise au Bureau international. 


Ce raDDort de recherche internationale comprend 3 feuilles. 




|~X~l II est aussi accompagne d une copie de chaque document relatif a I'etat de la technique qui y est cite. 


1. 


| | II a ete estim£ que certain es revendlcatlons nepouvalent pas falre I'objet d une recherche(voir le cadre I). 


2. 


[ | II y a absence d'unite de l'inventlon(voir le cadre II). 


3. 


| | La demande internationale contient la divulgation d'un listage de sequence de nucleotides oud'acldes amines et la 

recherche internationale a ete effectuee sur la base du listage de sequence 




□ 


depose avec la demande internationale 




□ 


fourni par le deposant separement de la demande internationale 

| | sans etre accompagnee d'une declaration selon laquelle il n'inctut pas d'elements 
allant au-dela de la divulgation faite dans la demande internationale telle 
qu'elle a ete deposee. 




□ 


transcrit par I'administration 


4. 


En ce.qui concerne le titre, [^] 


le texte est approuve tel qu'il a ete remise parle deposant. 




□ 


Le texte a ete etabli par I'administration et a la teneur suivante: 


5. 


En ce qui concerne I'abrege, 








le texte est approuve tel qu'il a ete remis parle deposant 




□ 


le texte (reproduit dans le cadre III) a ete etabli par I'administration conformement a la 
regie 38.2b). Le deposant peut presenter des observations a I'administration dans un delai 
d'un mois a compter de ia date d'expedition du present rapport de recherche internationale. 


6. 


La figure des desslns a publier avec I'abrege est la suivante: 




Figure n° 1 \%\ 
□ 


suggeree par le deposant. Q Aucune des figures 

... A , n'est a publier. 
parce que le deposant n a pas suggere de figure. K 




□ 


parce que cette figure caracterise mieux I'invention. 
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A. CLASSEMENT DE L'OBJET DE LA DEMANDE 

CIB 6 H04B7/005 H04L25/03 



Selon la classification intemalionale des brevets (C IB) ou 6 la tors selon la classification natiorale el la CIB 



B. DOMAINES SUR LESQUELS LA RECHERCHE A PORTE 



Documentation minim ale consultee (systeme de classification suivi des sym boles de classement) 

CIB 6 H04B H04L 



Documentation consultee autre que la documentation minimale dans la mesureou ces documents relevent des domaines surlesquels a porte la recherche 



Base de donnees electronique consultee au cours de la recherche internal ionale (nom de ta base de donnees, et si celaest realisable, termes de recherche 
utilises) 



C. DOCUMENTS CONSIDERES COMME PERTINENTS 



Categorie 0 Identification des documents cites, avec.le cas echeant, I'indicationdes passages pertinents 



no. des revendications visees 



EP 0 669 729 A (THOMSON CSF) 30 aout 1995 



ligne 1 
1 igne 20 



* abrege * 

voir page 2, ligne 1 - ligne 14 

voir page 3, ligne 58 - page 5, 

voir page 8, ligne 53 - page 9, 

voir revendication 1; figure 1 

US 5 479 446 A (M0UR0T CHRISTOPHE) 
26 decembre 1995 

* abrege * 

voir colonne 2, ligne 20 - ligne 39 
voir colonne 2, ligne 56 - colonne 4, 
1 igne 32 

voir figure UNIQUE 

-/— 



| )( j Voir la suite du cadre C pour la finde la liste des documents 



10 



Les documents de families de brevets sont indiquesen annexe 



0 Categories sped a les de documents cites: 

"A" document definissant I'etat general de latechnique, non 
considers com me particulierement pertinent 

"E" document anterieur, mais publie a la date dedepdt international 
ou apres cette date 

"L" document pouvant jeter un doute sur une revendcation de 
priorite ou cite pour determiner la date de publication d'une 
autre citation ou pour une raison speciale (telle qu'indiquee) 

"O" document se referant a une divulgation orale, a un usage, a 
une exposition ou tous autres moyens 

"P" document publie avant la date de depdtinternational, mais 
posterieurement a la date de priorite revendiquee 



T" document ulterieur publie apres la date de dep&t international ou la 
date de priorite et n'appartenenant pas a I'etat de la 
technique pertinent, mais cite pour com prendre le principe 
ou la theorie const it uant la base del'invention 

"X" document particulierement pertinent; ('invention revendiquee ne peut 
etre consideree comme nouvelle ou comme impliquant une activite 
inventive par rapport au document considere isolement 

"Y" document particulierement pertinent; I'invention revendiquee 

ne peut etre consideree comme impliquant uneactivite inventive 
lorsque le document estassocie a un ou plusieurs autres 
documents de meme nature, cette combinaison etant evidente 
pour une person ne du metier 

document qui fait partie de la meme famillede brevets 



Date a laquelle la recherche internationale a eteeffectivement achevee 



18 septembre 1998 



Date d'expedition du present rapport de recherche internationale 



29/09/1998 



Nom et adresse postal e de Tadministrationchargee de la recherche internationale 
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Fax: (+31 -70) 340-3016 



Fonctionnaire autorise 



Lopez Marquez, T 
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C.(sulte) DOCUMENTS CONSIDERES COMME PERTINENTS 



Catdgorie 0 Identification des documents cites, avec.le cas dcheant, I'indlcationdes passages pertinents 



no. des revendications visees 



EP 0 564 849 A (IBM) 13 octobre 1993 
* abrege * 

voir page 6, ligne 19 - ligne 48 
voir revendi cation 1; figures 1,2 
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FI 


934359 
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5432816 


A 
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Applicant's or agent's file reference 
FO 1002 WO 


See Notification of Transmittal of International 
FOR FURTHER ACTION p reIiminary Examination Report (Form PCT/IPEAM 1 6) 


International application No. 


International filing date {day/month/year) 


Priority date {day/month/year) 


PCT/FR98/00735 


10 April 1998(10.04.1998) 


14 April 1997(14.04.1997) 


International Patent Classification (IPC) or national classification and IPC 




H04B 7/005 






Applicant 


NORTEL MATRA CELLULAR 





This international preliminary examination report has been prepared by this International Preliminary Examining 
Authority and is transmitted to the applicant according to Article 36. 



This REPORT consists of a total of 



. sheets, including this cover sheet. 



□ This report is also accompanied by ANNEXES, i.e., sheets of the description, claims and/or drawings which have 
been amended and are the basis for this report and/or sheets containing rectifications made before this Authority 
(see Rule 70.16 and Section 607 of the Administrative Instructions under the PCT). 

These annexes consist of a total of sheets. 



This report contains indications relating to the following items: 
Basis of the report 
Priority 

Non-establishment of opinion with regard to novelty, inventive step and industrial applicability 
Lack of unity of invention 

Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

Certain documents cited 
Certain defects in the international application 
Certain observations on the international application 



I 




II 


□ 


III 


□ 


IV 


□ 


V 


El 


VI 


□ 


VII 




VIII 





Date of submission of the demand 

1 3 November 1 998 ( 1 3 . 1 1 . 1 998) 


Date of completion of this report 

09 July 1999 (09.07.1999) 


Name and mailing address of the IPEA/EP 
European Patent Office 
D-80298 Munich, Germany 

Facsimile No. 49-89-2399-4465 


Authorized officer 
Telephone No. 49-89-2399-0 
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al application No. 

PCT/FR98/00735 



I. Basis of the report 



1 . This report has been drawn on the basis of (Replacement sheets which have been furnished to the receiving Office in response to an invitation 
under Article 14 are referred to in this report as "originally filed" and are not annexed to the report since they do not contain amendments.). 

| [ the international application as originally filed. 

[^| the description, pages , as originally filed, 

pages , filed with the demand, 

pages , filed with the letter of 

pages , filed with the letter of 



^ the claims, 



Nos. 
Nos. 
Nos. 
Nos. 
Nos. 



1-14 



, as originally filed, 

, as amended under Article 19, 

, filed with the demand, 

, filed with the letter of 

, filed with the letter of 



the drawings, sheets/fig 
sheets/fig 
sheets/fig 
sheets/fig 



1/3-3/3 



, as originally filed, 
, filed with the demand, 
, filed with the letter of 
, filed with the letter of 



2. The amendments have resulted in the cancellation of: 

1 I the description, pages 

1 1 the claims, Nos. 



I 1 the drawings, sheets/fig 



3 I I This report has been established as if (some of) the amendments had not been made, since they have been considered 
— to go beyond the disclosure as filed, as indicated in the Supplemental Box (Rule 70.2(c)). 



4. Additional observations, if necessary: 
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V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 



Statement 
Novelty (N) 

Inventive step (IS) 
Industrial applicability (IA) 



Claims 
Claims 

Claims 
Claims 

Claims 
Claims 



2-14 



4-6, 9-14 



2-3, 7-8 



1-14 



YES 
NO 
YES 
NO 

YES 
NO 



Citations and explanations 

1. The following documents are referred to herein: 

Dl: EP 0 669 729 A 

D2: EP 0 665 665 A 

D3: EP 0 564 849 A 

Document D2 is not cited in the search report. 

2. Dl discloses (page 3, line 58 to page 4, line 17) a 
method wherein a multichannel transmission path (in 
the terms of the present application; in Dl, a 
channel consists of a plurality of paths) is 
estimated. The method requires the impulse response 
of the channels to be estimated and includes the 

.following steps: 

- acquiring a spatial statistic of the transmission 
path (page 4, line 8; page 5, lines 21-26), and 

- establishing a corrected impulse response by 
weighting said impulse response estimates with said 
spatial statistic and an estimate of the additive 
noise of said channels (page 4, line 8) . 
Therefore, the subject matter of claim 1 is not 
novel (PCT Article 33(2)). 
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The use, as the spatial statistic, of the 
correlation between pairs of transmission channels 
as well as the corresponding correlation matrix is 
part of the standard practice (cf. D2, formula 7, at 
the bottom of page 4). 

Therefore, the subject matter of claims 2 and 3 is 

not inventive (PCT Article 33(3)). 



4. D3 describes (figure 1; page 6, line 39; page 10, 
formula 14 and lines 14-16) a method for 
transmission channel estimation on the basis of a 
signal received by the channel, which received 
signal corresponds to a transmitted training 
sequence, characterised in that it comprises the 
steps of: 

- acquiring a statistic of the transmission channel, 
and 

- establishing a channel impulse response estimate 
weighted with said channel statistic by means of the 
received signal. 

Therefore, the subject matter of claim 7, in so far 
as it is dependent on claim 1 or 2, is not inventive 
(PCT Article 33 (3) ) . 

5. D3 describes (figure 2; page 6, line 39; page 10, 
formula 14; page 11, lines 18-24; page 14, line 29) 
a method wherein the statistic corresponds to an 
estimate of the impulse response covariance (matrix 
P) - 

Therefore, the subject matter of claim 8 is not 
inventive (PCT Article 33(3)). 

6. The subject matter of claims 4, 6 and 9 relates to a 
transmission path estimation method characterised in 
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that it comprises: 

(a) an estimation of the additive channel noise 
determined by normalising the energy of vector S- 
ACi, where S represents the received symbols 
corresponding to the training sequence, A is a 
measurement matrix constructed from said training 
sequence, and Ci is an estimate of the impulse 
response of the first transmission channel of the 
transmission path; 

(b) the corrected impulse response calculation 
formula indicated in claim 6, page 20, lines 6-9; 

(c) the weighted estimate calculation formula 
indicated in claim 9, page 21, lines 5-9. 

The prior art does not suggest the use of any of 
these three methods. Therefore, the subject matter 
of claims 4, 6 and 9 is novel and inventive. 

7. Claim 5 and claims 10-14 are dependent on claim 4 

and claim 9, respectively, and thus also comply, as 
such, with the requirements of novelty and inventive 
step of the PCT. 
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VII. Certain defects in the international application 



The following defects in the form or contents of the international application have been noted: 



2. 



It appears that the following corrections should be 
made to the original French text: 
page 4, line 7: " symboles" 
page 15, line 15: "a priori" 
page 20, line 26: "caracterise" 

Contrary to the requirement of PCT Rule 5.1(a) (ii) , 
the description does not indicate the relevant 
prior art disclosed in document Dl, and does not 
cite this document. 

The claims have not been drafted in two parts in 
accordance with PCT Rule 6.3(b). 
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VIII. Certain observations on the international application 

The following observations on the clarity of the claims, description, and drawings or on the question whether the claims are fully 
supported by the description, are made: 

As the reference signs placed in the claims are intended 
solely to facilitate the understanding thereof (PCT Rule 
6.2(b)), they cannot be used to restrict the scope of the 
claims. Therefore, the symbols Ci, C 2 , Ci, Cj, C n , N, G, G' , 
Nof X p/ A/i, X', v'i and S used in the formulae of the 
claims should have been defined in the claims. 
Furthermore, it should have been indicated that 1 
represents the identity matrix and that h denotes the 
Hermitian transform operator on the matrices. 
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(article 36 et regie 70 du PCT) 



PCT 



Reference du dossier du deposant ou du 
mandataire 

FO1002WO 


voir la notification de transmission du rapport d'examen 
POUR SUITE A DONNER preliminaire international (formulaire PCT/IPEA/416) 


Demande Internationale n" 
PCT/FR98/00735 


Date du depot international (jour/mois/annee) 
10/04/1998 


Date de priori te (jour/mois/annee) 
14/04/1997 


Classification internationale des brevets (CIB) ou a la fois classification nationale et CIB 
H04B7/005 


Deposant 

NORTEL MATRA CELLULAR et al. 



1 . Le present rapport d'examen preliminaire international, etabli par I'administaration chargee de I'examen preliminaire 
international, est transmis au deposant conformement a Tarticte 36. 

2. Ce RAPPORT comprend 6 feuilles, y compris la presente feuille de couverture. 

□ II est accompagne d'ANNEXES, c'est-a-dire de feuilles de la description, des revendications ou des dessins qui ont 
ete modifiees et qui servent de base au present rapport ou de feuilles contenant des rectifications faites aupres de 
('administration chargee de I'examen preliminaire international (voir la regie 70.16 et I'instruction 607 des Instructions 
administratives du PCT). 

Ces annexes comprennent feuilles. 



3. Le present rapport contient des indications relatives aux points suivants: 
I S Base du rapport 



II 


□ 


Priorite 


III 


□ 


Absence de formulation d'opinion quant a la nouveaute, I'activite inventive et la possibilite 
d'apptication industrielle 


IV 


□ 


Absence, d'un it e de I'invention 


V 


IS 


Declaration motivee selon Particle 35(2) quant a la nouveaute, I'activite inventive et la possibilite 
d'application industrielle; citations et explications a I'appui de cette declaration 


VI 


□ 


Certains documents cites 


VII 


ei 


Irregularis dans la demande internationale 


VIII 




Observations relatives a la demande internationale 



Date de presentation de la demande d'examen preliminaire 
internationale 

13/11/1998 


Date d'achevement du present rapport 

0 9. 07. 99 


Norn et adresse postale de I'admintstration chargee de 
I'examen preliminaire international: 

— — office europeen des brevets 

D-80298 Munich 

Tel. (+49-89) 2399-0 Tx: 523656 epmu d 
Fax: (+49-89) 2399-4465 


Fonctionnaire autorise 

Bauer. F vl ^ /) 

N" de telephone (+49-89) 2399 2737 ^U^^ 
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Demande internationale n> PCT/FR98/00735 



I. Base du rapport 

1 . Ce rapport a ete redige sur la base des elements ci-apres {les feuilles de remplacement qui ont ete remises a 
I'office recepteur en reponse a une invitation faite conformement a I article 14 sont considerees. dans ie present 
rapport, comme "initialement deposees" et ne sont pas jointes en annexe au rapport puisqu'elles ne contiennent 
pas de modifications.) : 

Description, pages: 

1 -18 version initiate 

Revendications, N": 

1 -14 version initiale 

Dessins, feuilles: 

1/3-3/3 version initiale 



2. Les modifications ont entraine I'annulation : 

□ de la description, pages : 

□ des revendications, n os : 

□ des dessins, feuilles : 

3. □ Le present rapport a ete formule abstraction faite (de certaines) des modifications, qui ont ete considerees 

comme allant au-dela de I'expose de I'invention tel qu'il a ete depose, comme il est indique ci-apres 
(regie 70.2(c)) : 

4. Observations complementaires. le cas echeant : 
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V. Declaration motivee selon l article 35(2) quant a la nouveaute, I'activite inventive et la possibility 
d'application industrielle: citations et explications a I'appui de cette declaration 

1. Declaration 



Nouveaute Oui : Revendications 2-14 

Non : Revendications 1 

Activite inventive Oui: Revendications 4-6,9-14 

Non : Revendications 2-3.7-8 

Possibility d'application industrielle Oui : Revendications 1-14 

Non : Revendications 



2. Citations et explications 
voir feuille separee 



VII. Irregularites dans la demande intemationale 

Les irregularites suivantes, concemant la forme ou le contenu de la demande intemationale, ont ete constatees : 
voir feuille separee 



VIII. Observations relatives a la demande intemationale 

Les observations suivantes sont faites au sujet de la clarte des revendications, de la description et des dessins 
et de la question de savoir si les revendications se fondent entierement sur la description : 

voir feuille separee 
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Concernant le point V 

Declaration motivee selon la regie 66.2(a)(ii) quant a la nouveaute, I'activite 
inventive et la possibility d'application industrielle; citations et explications a 
I'appui de cette declaration 

1 . II est fait reference aux documents suivants: 
D1: EP 0 669 729 A 
D2: EP 0 665 665 A 
D3: EP 0 564 849 A 

Le document D2 n'a pas ete cite dans le rapport de recherche. 



2. D1, page 3 ligne 58 - page 4 ligne 17, expose un procede incluant le sondage 
d'une voie de transmission formee de plusieurs canaux (au sens de la presente 
application. En effet, dans D1, un canal est compose de plusieurs voies), ce 
procede necessitant une estimation de la reponse impulsionnelle de ces canaux, 
et comprenant les etapes suivantes : 

- acquisition d'une statistique spatiale de cette voie de transmission (page 4 
ligne 8, page 5 lignes 21-26) 

- etablissement d'une reponse impulsionnelle corrigee par ponderation 
desdites estimations de reponses impulsionnelles au moyen de ladite 
statistique spatiale et d'une estimation du bruit additif desdits canaux (page 4 
ligne 8). 

L'objet de la revendication 1 n'est done pas nouveau (Article 33(2) PCT). 

3. ^utilisation, comme statistique spatiale, de la correlation des canaux de 
transmission pris deux a deux, ainsi que de la matrice de correlation 
correspondante, est une techique courante (cf D2, formule numero (7), au bas de 
la page 4). 

L'objet des revendications 2 et 3 n'est done pas inventif (Article 33(3) PCT). 

4. D3 (figure 1; page 6, ligne 39; page 10, formule 14 et lignes 14-16) decrit un 
procede de sondage d'un canal de transmission a partir d'un signal regu par ce 
canal, ce signal regu correspondant a une sequence d'apprentissage emise, 
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caracterise en ce qu'il comprend les etapes suivantes : 

- acquisition d'une statistique de ce canal de transmission 

- etablissement d'une estimation de la reponse impulsionnelle du dit canal 
ponderee par ladite statistique du canal au moyen du dit signal re<?u. 

L'objet de la revendication 7 comme dependant des revendications 1 ou 2 n'est 
done pas inventif (Article 33(3) PCT). 

5. D3 (figure 2; page 6, ligne 39; page 1 0, formule 1 4; page 1 1 , lignes 1 8-24; page 
14, ligne 29) montre un procede, dans lequel la statistique correspond a une 
estimation de la covariance (matrice P) de la reponse impulsionnelle. 

L'objet de la revendication 8 n'est done pas inventif (Article 33(3) PCT). 

6. L'objet des revendications 4, 6 et 9 porte sur un procede de sondage d'une voie 
de transmission caracterise respectivement par: 

a) une estimation du bruit additif du canal obtenue par normalisation de 
I'energie du vecteur (S-AC^, S representant les symboles re?us 
correspondant a la sequence d'apprentissage, A etant une matrice de 
mesure construite a partir de cette meme sequence d'apprentissage et 
une estimation de la reponse impulsionnelle du premier canal de 
transmission de la voie de transmission; 

b) la formule de calcul des reponses impulsionnelles corrigees indiquee a la 
revendication 6, page 20, lignes 6-9. 

c) la formule de calcul de I'estimation ponderee donnee dans la 
revendication 9, page 21, lignes 5-9. 

L'etat de la technique ne suggerait d'utiliser aucun de ces trois procedes. L'objet 
des revendications 4, 6 et 9 est done nouveau et inventif. 

7. Les revendications 5, respectivement 10-14 dependent des revendications 4, 
respectivement 9 et satisfont done egalement, en tant que telles, aux conditions 
requises par le PCT en ce qui concerne la nouveaute et I'activite inventive. 

Concernant le point VII 

Irregularites dans la demande intemationale 
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1 . II semble qu'il faille lire: 

page 4, ligne 7 "symboles" 

page 15, ligne 15 : "a priori" 

page 20, ligne 26 : "caracterise". 

2. Contrairement a ce qu'exige la regie 5.1 a) ii) PCT, la description n'indique pas 
I'etat de la technique anterieure pertinent expose dans le document D1 et ne cite 
pas ce document. 

3. Les revendications ne sont pas redigees en deux parties comme prevu par la 
regie 6.3b) PCT 



Concernant le point VIII 

Observations relatives a la demande intemationale 

Les signes de reference places dans les revendications etant uniquement destines a 
en faciliter la comprehension (Regie 6.2b) PCT), ils ne peuvent etre utilises pour en 
limiter la portee. Ainsi, les symboles (C 1f C 2f C if C if C n , N, G, G\ N 0 X p , A' h X\ v' if S) des 
formules des revendications auraient-ils du etre definis dans les revendications. 
D'autre part, il aurait fallu indiquer que I represente la matrice identite et que h designe 
I'operateur de transformation hermitienne sur les matrices. 
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(57) Abstract 

The invention concerns a method for estimating a 
communication path formed by several channels. Said 
method requiring an estimation of the impulse response 
function (Ci, C2, Co) of said channels comprises 
the following steps: acquiring a spatial statistic of said 
communication path; establishing a corrected impulse 
response function (d, C'2. — C'a) at least by weighting 
the pulse response functions (Ci, C% Cn) using the 
spatial statistic and an additive noise estimation (N01, N02. 
Nod) of these channels. 

(57) Abrege* 

L'invention conceme un proc6d6 de sondage d'une 
voie de transmission formee de plusieurs canaux. Ce 
proceae* qui necessite une estimation de la reponse im- 
pulsionnene (Ci, C2* Cn) de ces canaux comprend les 
etapes suivantes: acquisition d'une statistique spatiale de 
cette voie de transmission, etablissement d'une rfponse 
impulsionnelle corrigee (C*i, C'2. C'n) au moins par 
ponderation des estimations de ieponses impulsionneUes 
(Ci, C2. — , Cn) au moyen de la statistique spatiale et 
d'une estimation du bruit addiuf (Not, N02. -~» Nqq) de 
ces canaux* 
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Sondaae a ponderation spatiale d'une voie de transmission 
La presente invention concerne un procede de sondage 
d'une voie de transmission formee de plusieurs canaux. II 
s'agit done d'une technique connue sous le nom de diversite 
5 de reception , technique selon laquelle un recepteur comporte 
plusieurs antennes chacune associee & un canal de 
transmission distinct . En d 1 autres termes , 1 • invention 
propose une methode d 1 estimation de la reponse 
impulsionnelle de ces canaux de transmission. 

10 En effet, dans un systeme de transmission, notamment 

par ondes radio, le recepteur reyoit pour chaque canal de 
transmission un signal emis par un emetteur. Le signal emis 
subit des fluctuations d' amplitude et de phase dans le canal 
de transmission, si bien que le signal regu par le recepteur 

15 ne lui est pas identique. Les fluctuations du signal sont 
essentiellement dues a ce que l^omme de metier appelle 
1 1 interference intersymbole. Cette interference peut 
provenir de la loi de modulation employee pour la 
transmission et elle est egalement due a la propagation 

20 multi-trajets dans le canal. 

II se trouve que le signal reyu est generalement issu 
d'un grand nombre de reflexions dans le canal, les 
differents trajets empruntes par le signal emis coricfuisant 
ainsi a des retards varies au niveau du recepteur. La 

25 reponse impulsionnelle du canal repr6sente 1' ensemble de ces 
fluctuations, auxquelles est soumis le signal emis. II 
s f agit done la de la caracteristique fondamentale 
representant les transmissions entre l 1 emetteur et le 
recepteur . 

30 La reponse impulsionnelle du canal est utilisee 

notamment par un egaliseur qui a precisement pour fonction 
de corriger 1 1 interference intersymbole dans le recepteur. 
Une methode classique pour realiser une estimation de cette 
reponse impulsionnelle consiste a disposer dans le signal 

35 emis une sequence d 1 apprentissage formee de symboles connus. 
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Cette sequence est choisie en fonction de la loi de 
modulation et de la dispersion du canal, dispersion devant 
s« entendre ici comme le retard d'un symbole emis empruntant 
le trajet le plus long du canal par rapport a ce meme 
5 symbole empruntant le trajet le plus court. La dispersion 
est couramment exprimee comme un multiple de la duree qui 
separe deux symboles emis successifs, soit un nombre de 
"duree symbole" . 

A titre d'exemple, on citera deux techniques connues 
10 d' estimation de la reponse impulsionnelle d'un canal de 
transmission . 

La premiere technique fait appel a des sequences 
d'apprentissage particulieres dites sequences CAZAC, pour 
1* expression anglo-saxonne "Constant Amplitude Zero 

15 Autocorrelation". De telles sequences sont decrites dans 
!• article de A* MILEWSKI : "Periodic sequences with optimal 
properties for channel estimation and fast start-up 
equalization", IBM Journal of Research and Development, 
vol.27, N°5, Sept. 83, pages 426-431. 

20 systeme de radiocommunication cellulaire num§rique 

GSM fait appel a des sequences d'apprentissage TS formees de 
26 symboles notes a D a a 2 s prenant la valeur +1 ou -i. ces 
sequences possfedent les proprietes suivantes : 

20 

25 Yi a ± 2 - 16 

i = 5 

20 

2 a i a i + k = ° Si 0 < |k| ^ 5 

i = 5 

En notant d la dispersion du canal qui vaut 4 dans le 
cas du GSM, l 1 estimation de la reponse impulsionnelle prerid 
30 la forme d'un vecteur X a 5 composants notes x Q a x 4 . 

La sequence de symboles regus S correspondant a la 
sequence d ' apprentissage TS est formee elle aussi de 26 
symboles notes s Q a s 2 5- On suppose naturellement ici que 
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l'emetteur et le recepteur sont parf aitement synchronises et 

dans ce cas 1 1 estimation de la reponse impulsionnelle X est 

donnee par 1* expression suivante : 

20 

xk = — V ai s. pour 0 < k < 4 

16 i = 5 

5 La technique CAZAC presente l'avantage d'une grande 

simplicity de mise en oeuvre. Cependant, on remarque que 
chaque composante de la r6ponse impulsionnelle est etablie a 
partir de seulement 16 symboles regus. Etant donne que la 
sequence d 1 apprentissage comprend 26 symboles et que la 
10 dispersion du canal vaut 4, il y a de 1 1 information dans le 
signal regu qui n'est pas prise en compte et cela conduit a 
une reduction des performances par rapport a 1 1 ideal 
theorique . 

La deuxieme technique connue fait appel au critere des 
15 moindres carr6s. Elle est decrite notamment dans les 
demandes de brevet FR 2 696 604 et EP 0 564 849. En matiere 
de rappel, cette technique fait appel & une matrice de 
mesure A construite a partir de la sequence d 1 apprentissage 
TS de longueur n. Cette matrice comprend (n-d) lignes et 
20 (d+1) colonnes, d representant toujours la dispersion du 
canal. L'616ment figurant a la ieme ligne et a la jieme 
colonne est le (d+i-j)ieme symbole de la sequence 
d • apprentissage : 



25 



A = 



r a4 a3 a2 ai ao 
a 5 ^4 ^3 a2 ai 
as as a4 a3 a2 
*7 



^25 



a21 J 
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La sequence d 1 apprentissage est choisie telle que la 
matrice A fc A soit inversible ou l'operateur .t represente la 
transposition. C'est naturellement le cas pour les sequences 
CAZAC mais c'est egalement le cas pour d'autres sequences. 
5 Dans la sequence de symboles regus, on ne prend pas en 

compte les quatre premiers sq a S3 car ceux-ci dependent 
egalement de syrabole inconnus emis avant la sequence 
d 1 apprentissage, etant donne que la dispersion du canal vaut 
4 . Par un abus de langage on def inira done dorenavant le 

10 signal regu comme un vecteur S ayant pour composantes les 
symboles regus, S4, S5, S5, ... , S25. 

Des lors, 1' estimation de la reponse impulsionnelle 
prend la forme suivante : 

X = (At A)-l At . S 

15 Cette technique des moindres car res est un peu plus 

complexe que la precedente mais il faut noter que la matrice 
(A fc A)" 1 A fc est calcul6e une seule fois. On remarque ici que 
chacune des composantes de 1 1 estimation de la reponse 
impulsionnelle X est 6tablie a partir de 22 symboles regus 

20 et non pas de 16 comme dans le cas de la technique CAZAC. On 
doit done s'attendre a une amelioration des performances. 

Cependant, quelle que soit la technique utilisee, les 
reponses impulsionnelles de chacun des canaux de la voie de 
transmission sont considerees comme independantes . 

25 La pr§sente invention a ainsi pour premier objet un 

proc6de de sondage d'une voie de transmission qui prend en 
compte le fait que les differentes antennes sont 
spatialement li§es. 

Selon 1' invention, le procede de sondage d'une voie de 

30 transmission s 1 applique lorsque cette voie est formee de 
plusieurs canaux et il necessite une estimation de la 
reponse impulsionnelle C lf c 2/ .../ C n de ces canaux. Le 
procede comprend les etapes suivantes : 

- acquisition d'une statistique spatiale de cette voie de 
35 transmission, 
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- 6tablissement d'une reponse impulsionnelle corrigee au 
moins par ponderation de ces estimations de reponses 
impulsionnelles au moyen de cette statistique spatiale et 
d f une estimation du bruit additif de ces canaux. 
5 Avantageusement , cette statistique spatiale correspond 

a une estimation de la correlation des canaux de 
transmission pris deux a deux. 

Suivant un mode de realisation pr6f §rentiel , cette 
estimation de la correlation des canaux de transmission 

10 prend la forme d'une matrice de correlation spatiale ainsi 
definie que l 1 Element figurant a la ieme ligne et ct la jieme 
colonne soit obtenu par lissage du produit de la transposee 
hermitienne de la reponse impulsionnelle estimee du ieme 
canal et de la reponse impulsionnelle estimee du jieme 

15 canal. 

Selon une caractSristique additionnelle, un signal 
regu S par un premier canal correspondant & une sequence 
d ' apprentissage emise, en notant A la matrice de mesure 
associee a cette sequence d' apprentissage , 1' estimation du 
20 bruit additif Nqi de ce canal est obtenue par normalisation 
de 1 • energie du vecteur ( S - ACj ) • 

On peut prevoir que cette normalisation soit suivie 
d 1 une 6tape de moyennage . 

De plus, une matrice de bruit N 6tant formee a partir 
25 de l 1 estimation du bruit additif Ngi/ N 0 2/ N on des 

canaux, une matrice de ponderation spatiale G' etant ainsi 
definie a partir de la matrice de correlation spatiale G et 
de la matrice de bruit : 

G 1 = G (G + N)" 1 , 
30 les r§ponses impulsionnelles corrigees C'i, C*2/ ..♦ , C' n 
sont obtenues par 1 1 expression suivante : 
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Le proced6 de sondage d'une voie de transmission est 
done mis en oeuvre a partir de 1 'estimation de la r§ponse 
impulsionnelle des difffirents canaux consideres comme autant 
5 de canaux independants . Or les erreurs d' estimation sont 
inevitables. De maniere generale, la determination de la 
reponse impulsionnelle d'un canal de transmission unique est 
un probleme qui ne peut etre resolu de fagon exacte en 
presence de bruit additif. De plus les techniques 

10 anterieures font l'hypothese implicite que cette reponse 
impulsionnelle peut prendre une forme quelconque. 

La presente invention a ainsi pour deuxieme objet une 
mfethode d 1 estimation de la reponse impulsionnelle d'un canal 
de transmission qui presente une meilleure resistance au 

15 bruit additif ou, autrement dit, qui mene a une erreur 
reduite comparee a I 1 erreur d 1 estimation des techniques 
connues. Cette methode sera avantageusement appliquee pour 
1'un au moins des canaux formant la voie de transmission 
avant 1 • etablissement des reponses impulsionnelles corrigees 

20 par ponderation au moyen de la statistique spatiale et d'une 
estimation du bruit additif de ces canaux. 

Selon 1' invention, cette methode necessite un signal 
requ par un canal, ce signal regu correspondant a une 
sequence d'apprentissage emise. La methode comprend les 

25 stapes suivantes : 

- acquisition d'une statistique temporelle de ce canal de 
transmission , 

- etablissement de 1' estimation de la reponse impulsionnelle 
de ce canal, cette estimation etant ponderee par la 

30 statistique temporelle du canal au moyen du signal regu. 
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La statistique temporelle du canal represente une 
valeur de la r6ponse impulsionnelle anterieure a 
1 'acquisition du signal regu. La susdite ponderation 
introduit le fait que la response impulsionnelle afferente au 
5 signal regu a une valeur probablement plus proche de cette 
valeur anterieure qu'une valeur qui en serait tres eloignee. 
Ainsi, statistiquement , l'erreur d 1 estimation est diminuee. 

Avantageusement , cette statistique correspond a une 
estimation de la covariance de cette reponse impulsionnelle. 
10 Selon une premiere variante du proc^de, celui-ci 

comprend les etapes suivantes : 

- lissage de la reponse impulsionnelle et orthonormalisation 
au moyen d'une matrice de transformation W pour obtenir 
l 1 estimation de la covariance qui prend alors la forme d'une 

15 matrice L 1 , 

- recherche des vecteurs propres v^ 1 et valeurs propres X^ ' 
associees de cette matrice L 1 , 

- estimation de la reponse impulsionnelle instantanee du 
canal a partir du signal regu et application de cette 

20 matrice de transformation W pour former un vecteur X 1 , 

1 'etablissement de I 1 estimation ponderee Xp etant ainsi 
realisee : 



oil No est un nombre r6el strictement positif representant le 
25 bruit additif du canal. 

On peut ici prevoir que le bruit additif soit choisi 
egal & la plus petite des valeurs propres X± 1 . 

Lorsqu'un sous-ensemble de ces valeurs propres X±* 
presente une contribution inferieure a un seuil 
30 predetermine, on peut egalement prevoir que chacune de ces 
valeurs propres soit forcee a la valeur du bruit additif. 

La complexity s'en trouve reduite d'autant. 

Selon une seconde variante du procede, !• estimation de 
la covariance prenant la forme d'une matrice R, en notant A 
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la matrice de mesure associee a la sequence d 'apprentissage, 
1 'etablissement de 1 1 estimation ponderee est ainsi realist : 

Xp = (A fc A + NqR" 1 )" 1 At.S 
oil N 0 est un r6el strictement positif representant le bruit 
5 additif du canal. 

Par ailleurs, il est possible d'effectuer une etape 
d'orthonormalisation de la matrice R au moyen d'une matrice 
de transformation W pour obtenir une nouvelle matrice R 1 , 
1' estimation ponderee prenant alors la nouvelle forme 
1 0 sui vante : 

Xp = W fc (I + NqR 1 ~1 ) ~1 wt A !t .S 
ou la matrice A 1 est egale au produit de la matrice de 
transformation W et de la matrice de mesure A. 

Avantageusement , 1' expression (I + NqR'" 1 )" 1 est 
15 calculee au moyen du lemme d 1 inversion matricielle. 

La presente invention apparaitra maintenant de maniere 
plus detaillee dans le cadre de la description qui suit ou 
sont proposes des exemples de mise en oeuvre a titre 
illustratif, ceci en reference aux figures annexees qui 
20 represent ent : 

la figure 1, un diagramme identifiant les 
principales etapes d'un mode de mise en oeuvre du procede 
de sondage d'une voie de transmission selon 1' invention, 

- la figure 2, une premiere variante d' estimation de 
25 la reponse impulsionnelle d'un canal selon l f invention, 

- la figure 3, une deuxieme variante d' estimation de 
la r6ponse impulsionnelle d f un canal selon l 1 invention. 

L 1 invention sera presentee dans son application au GSM 
car ce systfeme a le merite djetre bien connu de I'homme du 
30 metier. II s'agit done la d'une presentation adoptee dans un 
souci de clarte mais il ne faut y voir en aucun cas une 
limitation de 1' invention a ce seul systeme. 

En reference a la figure 1, le procede de sondage 
d'une voie de transmission s' applique lorsque cette voie 
35 comporte au moins deux canaux de transmission, n canaux dans 
le cas general. Les canaux correspondent chacun a une 
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antenne distincte. On considere done que chaque canal a fait 
l'objet d'une estimation de sa reponse impulsionnelle 
respective Ci , C 2/ C n selon l'une quelconque des 

techniques disponibles . 
5 En premier lieu le procede prevoit 1 • acquisition d'une 

statistique spatiale de cette voie de transmission. Par 
statistique spatiale on entend un ensemble de donnees 
refletant le comportement de cette voie sur une periode 
pr^determinee que l'on conviendra d'appeler periode de 

10 correlation. En effet, les differentes antennes etant fixes, 
les signaux regus sur ces antennes presentent un certain 
degre de correlation. L 1 invention entend pr§cisement tirer 
parti de ce fait pour ameliorer la qualite d' estimation de 
la reponse impulsionnelle de l'un des canaux au moins. A 

15 titre d'exemple, une telle statistique peut etre obtenue au 
moyen d'une matrice de correlation spatiale G : 



20 



r v(c^c ± ) v(cjc 2 ) 



G = 



V( c^c n ) N 



VfC^) V(C^C 2 ) V < C 2 C n ) 



v(cj}c 1 ) v(cj}c 2 ) 



V < C n C n> 



oil 1 * operateur . h represente la transposition hermitienne . 

Cette matrice G carree de dimension (n,n) peut ainsi 
etre representee de maniSre generique par 1 1 Element 



9ij 



figurant a la ieme ligne et & la jieme colonne 
g. . = V(C^C. ) 



Cet element g^j est obtenu par lissage du produit cjrQ. 

25 au moyen des reponses impulsionnelles estimees C^, Cj des 
ifeme et jieme canaux obtenues durant la periode de 
correlation. Ce lissage est une estimation de la correlation 
de ces deux canaux. 
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On entend ici lissage dans un sens tres general, 
c ' est-a-dire toute operation permettant de lisser ou de 
moyenner le produit C^C-j sur la periode de correlation. 

Un premier exemple de lissage consiste a faire la 
5 moyenne de ce produit sur cette periode supposee comprendre 
e estimations successives de chacune des reponses 
impulsionnelles C^, C2/ * C n : 

v < c i c j> = ;S c h 

1 

10 Un deuxieme exemple de lissage consiste a actualiser, 

a la pieme estimation regue pour chacun des ieme et jieme 
canaux, I 1 expression de lissage _ ^C^Cj ) obtenue a la 

(p - l)ieme estimation au moyen d'un coefficient 
multiplicatif a, ce coefficient compris entre 0 et 1 etant 
15 souvent connu sous le nom de facteur d'oubli de lissage : 

V p (c i C j > = aC^Cj + (l-a)V p _ itcJCj ) 

L' initialisation peut se faire par tous moyens, 
notamment au moyen de la premiere estimation obtenue ou bien 
par une moyenne obtenue comme dans le premier exemple pour 

20 les premieres estimations regues. 

Le proced§ de sondage propose ensuite de realiser une 
estimation du bruit additif Nqi, Nq2/ •••/ Non present dans 
chacun des canaux au moyen des estimations des reponses 
impulsionnelles respectives C lf C 2 r •••/ C n de ces canaux. 

25 II sera propose differentes solutions pour mener a 

bien cette estimation du bruit et on traitera le cas d'un 
seul canal, le premier par exemple, §tant donne que les 
memes solutions s'appliquent a chacun des canaux. 

Une solution simple consiste a affecter Nqi d'une 

30 valeur pred6terminee qui reflete un seuil en dessous duquel 
il est peu probable que le bruit additif puisse descendre. 
Cette valeur pourrait etre determinee par une mesure de 
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rapport signal a bruit , ou par les performances du 
recepteur, ceci a titre d'exemple. 

Par ailleurs, 1 1 estimation du bruit Nqi peut etre 
obtenue au moyen de 1 1 estimation de la reponse 
5 impulsionnelle C± du premier canal et du signal regu 
correspondant S . Cette estimation du bruit a done pu etre 
realisee avant la mise en oeuvre du procede de sondage a 
ponderation spatiale selon quelque technique que ce soit. 
Toutefois, si ce n'etait pas le cas, on propose ici une 
10 methode appropriee lorsque l 1 estimation de la reponse 
impulsionnelle est acquise selon la technique des 

moindres carres . 

II s'agit done d'estimer directement le bruit additif 
a partir du signal regu S et de la ma trice de mesure A. En 
15 effet, en notant le vecteur bruit affectant le signal 

regu, il vient que : 

S = AC X + N x 

Compte tenu du fait que les vecteur s S et ont 22 
composanteS/ 1' estimation du bruit Nqi prend la forme 
20 suivante : 

N 0 1 = (^) (S » ACi)h (s - ACi) 

Naturellement cette estimation du bruit additif. Nqi 
peut etre moyennee ou lissee. 

A partir des estimations du bruit additif des 
25 differents canaux Nqi/ Nq2* N 0n/ on construit une 

matrice de bruit N dont 1* element figurant a la ieme ligne 
et & la ji&me colonne est not6 n^ j . 

Lorsque i est egal a j , cet 616ment n^ j prend la 
valeur du bruit estime Ngi« 
30 Lorsque i et j sont differents, 1 1 element n^j est nul . 

Par ailleurs, tou jours dans un souci de 
simplification, il est possible de calculer la valeur 
moyenne Nq des bruits estimes Nqi/ Nq2/ •*•/ Ng n e t de 
forcer chacun des elements njj de la diagonale de cette 
35 matrice N a cette valeur moyenne Nq. En notant I la matrice 
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identite, il s'ensuit que la matrice de bruit N prend la 
forme suivante : 

N = N 0 I 

La matrice de correlation G et la matrice de bruit N 
5 permettent de definir une nouvelle matrice, la matrice de 
ponderation spatiale G 1 : 

G 1 = G(G + N)" 1 
Les reponses impulsionnelles corrigees C^, C % 2 , •••/ 
C' n au moyen de cette matrice de ponderation spatiale G' 
10 sont ainsi definies : 















c 


















= G ' 


* 






• t 

n 


) 




• 

^ C n 


) 



ou l'op6rateur represente tou jours la transposition. 

L'une au moins des reponses impulsionnelles corrigees 
15 C'i est utilisee a la place de 1' estimation de la reponse 

impulsionnelle Ci dans le recepteur. 

A titre optionnel , 1 1 invention permet egalement 

d'ameliorer 1 1 estimation de la reponse impulsionnelle C l7 

^2/ •••/ C n de chacun des canaux, ceci avant la mise en 
20 oeuvre du procede de sondage a ponderation spatiale d'une 

voie de transmission, c'est-&-dire avant 1 ' etablissement de 

la ou des reponses impulsionnelles corrigees C' lf C f 2 / .../ 

C'n- 

On prevoit ainsi 1 1 acquisition d ' une statistique 
25 temporelle pour l'un au moins de ces canaux, le premier par 
exemple. Par statistique temporelle, on entend un ensemble 
de donnees refletant le comportement de ce canal considere 
independamment des autres canaux sur une periode 
predeterminee que l'on conviendra d'appeler periode 
30 d 1 analyse. II s'agit done d"une representation du 
comportement moyen du canal pendant la periode d' analyse. 
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Cette statistique peut etre etablie par quelque moyen que ce 
soit et en quelque lieu que ce soit. En effet, 
1 •etablissement de la statistique peut prendre place dans un 
§quipement quelconque du reseau de radiocommunication . Ce 
5 qui importe c'est que le recepteur puisse acquerir cette 
statistique* 

A titre d'exemple, une telle statistique temporelle 
peut etre obtenue de la maniere suivante. 

Au moyen du signal requ S pendant la p6riode 
10 d' analyse, on calcule une estimation X de la reponse 
impulsionnelle selon une methode connue. 

Si l'on retient la technique des moindres carres cette 
estimation X vaut : 

X = (A fc A)" 1 At.S (1) 
15 On rappelle ici que l'emetteur et le recepteur sont 

supposes synchronises a mieux que un demi-symbole pres, 
auquel cas le signal de reception est forme par le vecteur S 
dont les composantes sont les symboles regus S4 a S25 
synchrones des symboles a4 a a 2 s de la sequence 
20 d'apprentissage TS. Si une telle synchronisation n'etait pas 
acquise, plusieurs solutions sont disponibles pour 
1 ' acquerir et on en citera deux a titre d 1 exemple , 

La premiere solution consiste a decaler le signal regu 
en avance ou en retard de j periodes symboles de sorte que 
25 Sjt = (s 4 _j, ss.j, s 6 _j f s 2 s-j). 

On calcule alors 1 1 estimation Xj pour chacun des 
vecteurs Sj et on retient la valeur pour laquelle Xj h -Xj 
est maximal. Cette valeur jjj donne la synchronisation 
escomptee et il suffit de remplacer le vecteur S dans 
30 1' equation (1) par le vecteur Sj^* 

La deuxieme solution consiste a augmenter 
artif iciellement la dispersion d du canal d'une quantite "2q 
predeterminee . On peut alors definir une matrice de mesure 
modifiee A m comprenant (n-d-2q) lignes et (d+2q+l) colonnes, 
35 En reprenant les valeur s 26 et 4 respectivement pour n 
et d : 
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Am = 



+ 2g • • 
a 5 + 2q • • 
a 6 + 2q • • 
a 7 + 2q •• 



a4 a3 a 2 ai ao 
a 5 a 4 a3 a2 &i 
a6 as a4 a3 a2 



a 25 



a 21 - 2q 



II faut 

alors reduire le nombre de composantes du signal recu S de 
cette meme quantite 2q et on peut convenir de retenir le 
5 vecteur S m modifie : 

s m t = ( s '4+q/ s' 5+q/ ... s' 2 5-q) 
On obtient ainsi une estimation modif iee X m : 

x m = (Am* A m ) A m t.s m 
Cette estimation modifiee X m comporte d+2q+l 
10 composantes : 

Xm* = ( x -q/ . .. Xq, xi, • ••/ ... x 4+ g) 

On recherche alors la valeur j M de j comprise entre -q 

et +q qui maximise 1' expression suivante : 
4 

Z xj + k xj +k 

k = o 

15 ou l'operateur .* represente la conjugaison complexe. 

Cette valeur j M determine 1' estimation X de la reponse 
impulsionnelle pour une dispersion d=4 : 

xt = ( x jM' x jM+l/ x jM+4) 
La synchronisation s'en deduit immediatement en 
20 appliquant le decalage j M au signal regu S. 

On peut des lors appliquer a nouveau 1 • equation ( 1 ) . 
On construit maintenant une matrice de lissage L par 
lissage des differentes estimations X obtenues pendant la 
periode d' analyse pour obtenir une estimation de la 
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covariance associee a cette reponse impulsionnelle. On 
entend ici lissage dans un sens tres general, c'est-a-dire 
toute operation permettant de lisser ou de moyenner la 
reponse impulsionnelle sur la periode d' analyse. On obtient 
5 ainsi une representation statistique du comportement du 
canal de transmission. Ce lissage peut etre obtenu notamment 
par l'une des deux methodes proposees ci-dessus, 
1 ' expression la plus simple de cette matrice de lissage L 
§tant la suivante : 
io L = ^]Txx h 

ou m correspond au nombre de sequences d 1 apprentissage sur 
lequel est calcule le lissage. 

On admet ici que cette matrice de lissage peut etre 
approchee par 1 1 equation suivante : 

15 L ~ (AtA)" 1 N 0 + R (2) 

ou No represente tou jours le bruit present dans le 
canal de transmission ou bruit additif et ou R est une 
matrice que l'on a coutume d'appeler statistique a priori du 
canal car elle represente le comportement du canal 

20 abstraction faite du bruit. 

On admet egaleraent que la matrice de mesure A est bien 
conditionnee , c'est-a-dire que les valeurs propres de la 
matrice A^a sont tres proches les unes des autres. Dans ce 
cas, il est interessant de proc§der & 1 'orthonormalisation 

25 des vecteurs constitu6s par les lignes de la matrice de 
mesure A, mais il ne faut pas voir la une limitation de 
1 1 invention . 

Pour ce faire, on emploie une matrice de 
transformation W telle que : 
30 A' = AW et A'tA 1 = I 

oCi I represente la matrice identity. 

En notant L f la matrice ainsi definie : 
L = WL'wt, 

il vient que liquation (2) peut maintenant s f 6crire : 
35 L 1 - N 0 I + R f (3) 
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Selon une premiere variante representee a la figure 2, 
on remarque que les vecteurs propres vi* de L' et de R' 
sont identiques tandis que les valeurs propres de L' et 
X± R* sont decalees de N 0 . soit en prenant tou jours la 
5 meme valeur de 4 pour la dispersion du canal, pour tout i 
compris entre 0 et 4 : 
V ± ' = v ± 
Xi« = Xi + N 0 

II apparalt ainsi que la determination des vecteurs 

10 propres et valeurs propres de R 1 et celle de L' sont 

identiques sous reserve que Nq soit connu. 

L'etape d' estimation du bruit sera decrite plus loin 

pour rendre 1' expose plus clair, bien que cette etape 

precede celle qui va maintenant etre explicitee. 

15 II convient done maintenant de rechercher les couples 

valeur-propre/vecteur-propre pour l'une ou 1' autre des 

matrices L' ou R 1 . Cette etape ne sera pas plus detail lee 

car bien connue de l'homme du metier. Par ailleurs, il va 

sans dire que l'on peut annuler les valeurs propres dont la 

contribution est jugee non significative. Par exemple, si 

ces valeurs propres sont classees par ordre decroissant, on 

supprime les dernieres qui sont telles que leur somme soit 

inferieure a un seuil predetermine. 

L'etape suivante consiste a estimer la reponse 

25 impulsionnelle instantanee X etablie selon l'une quelconque 

des techniques connues a partir du signal recu correspondant 

a la derniere sequence d ' apprentissage regue. En notant 

X = WX', cette derniere estimation est pondered par la 

methode suivante pour obtenir une ponderation temporelle Xp 

30 de la reponse impulsionnelle instantanee : 

4 ( k ^ 
Xp = T i ( V *X») 

i = 0 V 1 0 j 

4 ( \ ' - M "\ 



20 



W V. 



i = 0 



W V. 

1 
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Pour obtenir la ponderation Xp il faut done estimer le 
bruit additif Nq» 

Cette estimation du bruit peut §tre effectuee 
notamment selon l'une des methodes exposees plus haut. 
5 Une autre possibility consiste a considerer que la 

derni&re valeur propre, (la plus faible) de la matrice de 
lissage L est egale a Nq . 

' = Nq ou ^4 = 0. 
Quelle que soit la m6thode retenue, la ponderation 
10 temporelle Xp de 1 1 estimation de la r6ponse impulsionnelle 
instantanee peut alors etre realis6e comme mentionne ci- 
dessus * 

Selon une seconde variante representee a la figure 3, 

on 6tablit 1 1 estimation ponder§e Xp directement comme suit : 

15 Xp = (A fc A + NQR-ij-iAt.S 

ou bien en reprenant la matrice de transformation W 

definie ci-dessus : 

Xp = W (I + NqR 1 " 1 )" 1 W fc A l1: .S (4) 

Conf ormement a 1 1 equation ( 3 ) : 

20 R' = L 1 - N 0 I 

II convient la aussi d 1 estimer le bruit additif Nq. 

Une solution avantageuse pour obtenir la ponderation 

temporelle Xp consiste ci adopter la m§thode suivante. 

On divise la matrice R * par Nq : 
R 1 

25 B = 

N 0 

II s'ensuit que : 

I + NqR'" 1 = I + B" 1 
On utilise le lemme d' inversion matricielle pour 
calculer la matrice de ponderation P = (I + B" 1 )" 1 . 
30 Ainsi, en notant e^ les vecteurs canoniques, on 

procede a 1 1 iteration suivante : 

- initialisation : 

P = B 

- pour i variant de 0 a d (4 dans le cas present) : 
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P = 



P - 



1 + e^Pe. 
1 l 



P 6tant connu, il reste a etablir la ponderation Xp 
selon 1 1 equation ( 4 ) . 

On remarquera que la matrice de ponderation P n f est 
5 pas necessairement calculee a I 1 occasion de chaque nouvelle 
sequence d ' apprentissage 6mise. Elle peut etre calculee a un 
rythme plus lent car elle varie sensiblement au meme rythme 
que R 1 et done plus lentement que le signal regu S. 



10 realisee sans avoir recours a la reponse impulsionnelle 
instantanee. Elle est produite directement a partir du 
signal regru S. 

Quelle que soit la variante retenue, e'est done 
1 * estimation ponderee Xp qui sera avantageusement utilisee 

15 coirane estimation de la reponse impulsionnelle pour la 

mise en oeuvre du precede de sondage a ponderation spatiale 
d'une voie de transmission, e'est-a-dire pour produire une 
ou plusieurs reponses impulsionnelles corrigees. 



On remarquera egalement que l 1 estimation ponderee est 
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REVINDICATIONS 

1) Procede de sondage d'une voie de transmission 
formee de plusieurs canaux, ce proc6de necessitant une 
estimation de la reponse impulsionnelle (C^, C2, /C n ) de 

5 ces canaux, caracterise en ce qu'il comprend les etapes 
suivantes : 

- acquisition d'une statistique spatiale de cette voie de 
transmission, 

- etablissement d 1 une reponse impulsionnelle corrigee ( C 1 1 , 
10 C r 2/ / c, n) au moins par ponderation desdites estimations 

de reponses impulsionnelles (C^, C2/ f C n ) au moyen de 

ladite statistique spatiale et d'une estimation du bruit 
additif (N 0 i/ N 0 2/ •••/ N 0n) desdits canaux. 

2) Procede selon la revendication 1, caracterise en ce 
15 que ladite statistique spatiale correspond a une estimation 

de la correlation desdits canaux de transmission pris deux a 
deux. 

3) Procede selon la revendication 2, caracterise en 
ce que cette estimation de la correlation des canaux de 

20 transmission prend la forme d'une matrice de correlation 
spatiale (G) ainsi definie que 1 'element (gij ) figurant a la 
ieme ligne et a la jifeme colonne soit obtenu par lissage du 
produit (C^Cj) de la transposee hermitienne de la reponse 

impulsionnelle estimee (C^) du i&me canal et de la reponse 
25 impulsionnelle estimee (Cj) du jieme canal. 

4) Proc6de selon la revendication 3 caracterise en 
ce que, un signal regu S par un canal correspondant ct une 
sequence d'apprentissage 6mise, en notant A la matrice de 
mesure associee a ladite sequence d 'apprentissage, 

30 1' estimation du bruit additif (Nqi) de ce canal est obtenue 
par normalisation de l'energie du vecteur (S - AC^). 

5) Proc6de selon la revendication 4, caracterise en ce 
que ladite normalisation est suivie d'une etape de 
moyennage . 

35 6) Procede selon I'une quelconque des 

revendications 3 a 5 caracterise en ce que, une matrice de 
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bruit (N) etant form6e a partir de 1' estimation du bruit 
additif (N 0 i/ N 0 2/ •••/ N 0n ) des canaux, une matrice de 
ponderation spatiale (C ) etant ainsi definie a partir de 
ladite matrice de correlation spatiale (G) et de ladite 
5 matrice de bruit : 

G' = G (G + N)" 1 , 

lesdites reponses impulsionnelles corrigees (C' lf C ! 2/ / 

C' n ) sont obtenues par 1' expression suivante : 
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10 7) Procede selon l f une quelconque des 

revendications 1 a 6 caracterise en ce que, avant 
1 f 6tablissement de ladite reponse impulsionnelle corrigee 
(C 1 !) d'un canal , le signal regu (S) par ce canal 
correspondant a une sequence d'apprentissage emise, il 

15 comprend les etapes suivantes : 

- acquisition d'une statistique temporelle de ce canal de 
transmission , 

- etablissement de 1» estimation de la reponse impulsionnelle 
(Ci) dudit canal, cette estimation (Xp) etant ponderee par 

20 la dite statistique temporelle du canal au moyen dudit 
signal regu ( S ) . 

8) Procede selon la revendication 7 , caracterise en 
ce que ladite statistique temporelle correspond a une 
estimation de la covariance de ladite reponse 

25 impulsionnelle. 

9) Procede selon la revendication 8 caracterise en 
ce qu'il comprend les etapes suivantes : 

lissage de ladite reponse impulsionnelle et 
orthonormalisation au moyen d'une matrice de transformation 
30 W pour obtenir ladite estimation de la covariance qui prend 
alors la forme d'une matrice L 1 , 
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- recherche des vecteurs propres (v^ 1 ) et valeurs propres 
(Xj/) associees de cette matrice L* , 

- estimation de la reponse impulsionnelle instantanee du 
canal a partir dudit signal regu (S) et application de cette 
matrice de transformation W pour former un vecteur X' , 

1 1 etablissement de ladite estimation ponderee (Xp) etant 
ainsi realisee : 



= z p_p (v . h xi) j Bv . h 



oti Nq est un nombre reel strictement positif representant le 
10 bruit additif dudit canal. 

10) Procede selon la revendication 9, caracterise en 
ce que ledit bruit additif (Nq) est choisi egal a la plus 
petite desdites valeurs propres (^i*). 

11) Procede selon l'une quelconque des revendications 9 
15 ou 10 , caracterise en ce qu*un sous-ensemble desdites 

valeurs propres (A^') presentant une contribution inferieure 
a un seuil predetermine, chacune de ces valeurs propres est 
forcee a la valeur dudit bruit additif (Nq). 

12) Procede selon la revendication 8 caracterise en 
20 ce que ladite estimation de la covariance prenant la forme 

d'une matrice R, en notant A la matrice de mesure associee a 
ladite sequence d'apprentissage, 1 1 etablissement de ladite 
estimation ponderee (Xp) est ainsi realist : 
Xp = (A*- A + NqR" 1 )" 1 A fc .S 
25 oil Nq est un reel strictement positif representant le bruit 
additif dudit canal . 

13) Procede selon la revendication 12, caracterise 
en ce qu'il comprend une §tape d ' orthonormalisation de 
ladite matrice R au moyen d'une matrice de transformation W 

30 pour obtenir une nouvelle matrice R' , cette estimation 
ponderee prenant alors la nouvelle forme suivante : 

Xp = wt(I + NqR 1 " 1 )- 1 W fc A ft .S 
od la matrice A' est 6gale au produit de cette matrice de 
transformation W et de ladite matrice de mesure A. 
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14) Procede selon la revendication 13, caracterise 
en ce que 1' expression (I+NqR'" 1 )" 1 est calculee au moyen du 
lemme d 1 inversion matricielle. 
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L'(XX n ) 



L , (XX n ) = N 0 l + R' 



Nq (estimation) 



(Vj.A'j)? i = 0,...,d 
Lv'j = X\ v'j 



X = (A t A)- 1 A t .S 



WX'=X 



Xp=S £l ? !D (v , [ 1 X t )Wv'i 
i=0 



FIGURE 2 
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L'(XX n ) 



L'(XX n )^NQl + R' 



N Q (estimation) 



P = (1 + NqR'" 1 )~ 1 



A'= AW 



X n = WPA'* S 



FIGURE 3 
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